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We investigate the origin of the transition inside the desynchronization state via 
phase jumps in coupled chaotic oscillators. We claim that the transition is governed 
by type-I intermittency in the presence of noise whose characteristic relation is (Z) oc 
exp(a|e( — e^/^) for — e < and (l) oc {et — e)^^/^ for — e > 0, where (/) is the 
average length of the phase locking state and e is the coupling strength. To justify 
our claim we obtain analytically the tangent point, the bifurcation point, and the 
return map which agree well with those of the numerical simulations. 
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I. INTRODUCTION 

Recently synchronization phenomena in coupled chaotic oscillators have been extensively 
studied because of their fundamental importance in areas of science and technology such 
as laser dynamics, electronic circuits, and biological systems|l|, Q, 0]. Due to interaction 
between two coupled chaotic oscillators, various features of synchronization are observed 
depending on the coupling strength. For example, when two identical chaotic oscillators 
are coupled, they can be synchronized either perfectly or intermittently jsl. In another 
case, where they are coupled with slight parameter mismatch, non-synchronization, phase 



synchronization j6|, Q, 0], or lag synchronization is observed depending on the couphng 
strength. Among these features, the noteworthy one of our study is phase synchronization 
(PS)p, 0, 0]. As is widely understood, above the critical strength of the coupling for the 
transition to PS, suitably defined phases of two chaotic oscillators are locked while their 
amplitudes remain chaotic and uncorrelated. And below the critical value, the phases of 
two oscillators are intermittently unlocked, that is, 27r phase jumps interrupt the phase 
locking states irregularly. So it can be said that the transition from nonsynchronous state 



to PS state is typically accompanied by an intermittent sequence of 2tt phase jumps 



To explain the transition mechanism to PS, two different approaches have been intro- 
d„ced: .opo>o,calBBB a.d .aUsHcalflQ. The t„po.og,ca. approach — tKat 

the behavior of coupled chaotic oscillators is analogous to that of a chaotic oscillator driven 



by an external chaotic signaljl^. From this study it was concluded that the phase jump 
phenomenon stems from a boundary crisis mediated by an unstable- unstable pair bi- 
furcation, which is termed eyelet intermittencv jl^ . Meanwhile, the statistical approach 
focuses on a phase equation which describes the phase difference of coupled chaotic oscilla- 
tors. A potential modulated by multiplicative noise was derived by the analysis of the phase 
equation and it was concluded that the transition mechanism is similar to that of eyelet 
intermittencv[lo| . 

By taking the statistical approach, however, we are to show in this paper that the charac- 
teristic relation of the average length of the phase locking state (or average laminar length) 
on the coupling strength before the transition to PS in coupled Rossler oscillators follows 
the relation of type-I intermittency in the presence of noise. To validate our argument, we 
derive a one-dimensional phase equation from the coupled Rossler oscillators and obtain 
the bifurcation point, the tangent point that is a 7i/2 phase locking state, and the return 
map analytically, which exactly agree with those of the numerical simulations. And then 
we obtain the characteristic relations of (/) ~ exp(Q;|et — e|^/^) for (cf — e) < and of 
(/) ~ {et - e)-i/2 for fe.- e) > which are the same as those of type-I intermittency in the 
presence of noise 13], where (/) is the average length of the phase locking state, e is the 
coupling strength and is the tangent bifurcation point. 

In our approach, the analytic scaling rule is obtained from the Fokker-Planck equation 



22j | and so it describes the asymptotic behavior near PS regime in which the effective noise 



can be reasonably approximated to be Gaussian 
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For this reason, it seems that 



3 



critical coupling does not appear in our formalism. However in real systems there is one 
critical coupling ec which is defined as a border of existence of phase slips, because of a finite 
amplitude of the effective noise. 

In section II, we will obtain phase difference equation of coupled Rossler oscillators from 
one-dimensional phase equation and compare the analytical and numerical return maps of 
the phase difference in Section III. After that we will obtain the characteristic relations 
of intermittent phase locking time according to the coupling strength by Fokker-Planck 
equation in Section IV, discuss the results in Section V, and conclude our study in Section 
VI. 



II. ANALYTIC STUDY OF COUPLED ROSSLER OSCILLATORS 

The transition to PS was first observed in coupled Rossler oscillators which have a slight 
parameter mismatch 

ii,2 = -^^1,22/1,2 - 2:1,2 + e(a;2,i - a;i,2), 

yi,2 = t^i,2a;i,2 + 0.15?/i,2, (1) 
ii,2 = 0.2 + zi.2(a;i,2- 10.0), 

where cui^ (= 1.0±0.015) are the overall frequency of each oscillator and e is the coupling 
strength. By transforming the above equation to the polar form, we obtain the following 
one-dimensional phase equation, which describes the phase difference of the two oscillators: 



where, 



= Aa; + A{ei, 62, e) sin0 + B{e,, 62)6 + ^(^i, ^2), (2) 



A(^i,^2,e) = (e + O.15)cos(^i + 02)- + 



2 Ri R' 



■2 



m,e2) = ^sin(^i)-^sin(^2), (3) 
rCi ri2 



Auj = uji —UJ2, = 6*1 — 62, 6*1,2 = arctan(yi, 27x1,2), and i?i,2 = \Jx\,2 + ?/i,2- Here 5, ^, and 
(e + 0.15) cos(6^i + O2) in A are fast fiuctuating terms in comparison with slowly varying 0. 
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FIG. 1: The values of {A), {B), (^), and e according to the coupling strength: (a) the mean 
values of A(circle), B(square), and ^(filled circle); (b) and (c) the values of ^ and e (the arrows 



point out the bifurcation point e^). 



Whereas, in the previous studies, the authors neglected all of the fast fluctuating 
terms^^, we found that they play a crucial role in the analysis of the transition mech- 
anism. Equation (2) can be transformed into the following simple form: 

^ = Acu+(A)sin0 + e, (4) 

where | = - (0) + {B - {B))e + {A - {A)) sin0, and AcD = Acj + + {B)e. Here {A) 
and {B) are the mean value of A and B respectively. This equation is similar to the one 
describing a phase locking of the periodic oscillator in the presence of noise [17j. 

In Eq. (4), if we turn off ^, the analysis of the system stability is straightforward. If 
^ = the system ends time-evolution and remains at 0*, where 

= arcsin(-^). (5) 

Here the condition of 0* being stable is |-^| < 1. Tangent bifurcation occurs at ^ = — 1 
and the tangent point is 0^ = ±| ± 2?™, where n = 0, 1, 2, ... and the sign ± depends on 



Auj. In our system since Au is positive, on 
near ±| ± 2nn and 27r phase jumps occur 



i*_ appears. This explains why is locked 



1( 
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III. NUMERICAL STUDY AND RETURN MAPS 

If Eq. (4) is expanded around the tangent point (t)\, the following equation is obtained: 
^ e + a02 where = 0- 7r/2, e = Au; + {A), and a = -{A)/2. Here if | is absent, 
moves very slowly around the tangent point (p*. So the dynamics of is mainly governed by 
^ in the situation |^| ^ |e|. Then we can regard 0^ as a constant when the two oscillators 
are in a locked state. We obtain a local Poincare map by integrating the above equation 
during the period that oscillator 1 completes every N rotation (the structure of the local 
Poincare map is invariant with respect to the number of N as far as N is small enough in 
comparison with the average length of the phase locking state). The local Poincare map is 
given by: 

(pn+l = 4>n + ^' + a(j)l + ^n, (6) 

where e' = 6(T„), a = Adt), and = JZ. idt + Adt - (/;"_^ Adt)). Here T„ 

1, where r„ is the overall time that oscillator 1 takes to complete N rotations. This 
is the very local Poincare map of type-I intermittency in the presence of noise when ^„ 
acts as random noise. In this equation tangent bifurcation occurs at the point e' = 0, which 
meets the condition of ^ = — 1. To determine the tangent bifurcation point, we calculate 
{A), {B), {^), and e according to the coupling strength e as presented in Fig. 1 (a), (b), 
and (c), respectively. The figures show that the value of e at the bifurcation point is 0.0277. 

In order to verify the above analytic results, we construct the return maps directly from 
Eqs. (l)-(3) by obtaining (p for = 1. Figure 2(a), (b), and (c) show the return maps 
before, near, and after the tangent bifurcation point, respectively. The figures show that as 
the coupling strength e increases, the shadow curve approaches the diagonal line. When the 
curve touches the diagonal line, tangent bifurcation occurs and the tangent point is | (see 
the inset arrow in Fig. 2(b)), which agrees well with the one obtained in the above. 

We define a measure A = | (n+ — n_)/ {n^ + n_)\, where n+ is the number of points above 
the diagonal line in the total number of points inside of the rectangle of x / and n_ is 
the number of points below the diagonal line. Then the measure shows the average ratio 
between above and below the passage near the tangent point. So the minimum value of A 
indicates the bifurcation point. Figure 2(d) shows a sharp minimum at e ^ 0.0277. This 
value again agrees with what we obtained from Fig. 1(b) and (c). 
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FIG. 2: The return maps {4>n,(pn+i) and A'^ as a function of e: (a), (b), and (c) are the return 
maps before (e = 0.018), near (e = 0.028), and after (e = 0.035) the tangent bifurcation point, 
respectively; (d) shows that the minimum value of A appears at et = 0.0277 (when k = 0.03 and 
/ = 0.6). 



The shadow curves are well fitted to the following form of type-I intermittencv|18l. I20| : 



'->n+l 



+ «0„ + e, 



(7) 



where 



2 ' 



0.094, and e = et — e. The coefficient a in Eq. (7) agrees well with 



a in Eq. (6), since (T„) ^ 2tt and the mean value of A is —0.03. (The overall frequency 
Ui^2 ~ 1-0.) This confirms that the phase equation of the coupled Rossler oscillators coincide 
with the structure of type-I intermittency. 



IV. RESULT FROM FOKKER-PLANCK EQUATION 

Under the long laminar length approximation, Eq. (6) can be transformed into the 
differential form ^ = e' + + ^. Then we can obtain the Fokker-Planck equation (FPE) 



by regarding ^„ as Gaussian white noise 



21 



The probability distribution and auto- 



correlation of C,n are examined for various A^s. So we find they are in better agreement 
with the Gaussian distribution and 5-correlation, respectively, as becomes larger. Figure 
3 shows the probability distribution and auto-correlation of ^„ for = 25 which coincide 
with the Gaussian profile with a dispersion of 0.4 and the 5-function, respectively. (A^ = 25 
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FIG. 3: The distribution and auto-correlation of for N = 25. (a) the distribution of well 
fitted to the Gaussian profile, and (b) the auto-correlation of ^„ almost 5-correlated. 

is still much less than the average length of the phase locking state that we have obtained.) 
The noise with the Gaussian distribution is not bounded whereas ^„ in coupled Rossler 
oscillators is bounded within (—1.3, 1.3). However, from the numerical data, we find that 
the probability of the occurrence of events outside the bounded region of ^„ is less than 



1.0 X 10 ^%. So it is a negligible effect on the average length o: 



the phase locking state. 



pn, 



From the FPE with appropriate boundary conditions [15|, ll6|, ll9|, |2l[, we can obtain the 
characteristic form of the average laminar length according to the coupling strength e as 
follows: 

(/) = (/o)exp(a|ei-e|3/2), (8) 

where a is the constant and (/q) is the average length of the phase locking state at the tangent 
bifurcation point. Figure 4(a) shows the average length of the phase locking state for = 1 
as a function of \et — e\ in the region 0.0200 < e < 0.0302. The slope in the space In \et — e\ 



versus ln(ln(/) — ln(/o)) is 1.44 as shown in Fig. 4(b) and its inset (c). The line fits wel 



within 



4.0% error from the 3/2 slope. The slope of the tail in Fig. 4 (b) converges to 1 25[ which 
shows the transient regime from (1) oc (e^ — e)~^/^ to (/) oc exp(a|et — ep^^). Figure 4(d) is 
the plot of ln(e( — e) versus ln(/) in the region 0.024 < e < 0.0277 and the tail also clearly 
shows the transient regime. The straight line fits well with the —1/2 slope. This means that 
the characteristic relation deforms from the conventional scaling rule (/) oc (e^ — e)^^/^ to 
(/) oc exp(a|et — e|^/^), as the coupling strength crosses the tangent bifurcation point. Thus 
we can understand that the average length of the phase locking state agrees well with the 
characteristic relation of type-I intermittency in the presence of noise not only in the region 
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ln\E-E\ ln\&-e\ 

FIG. 4: The numerical verification of the characteristic relation: (a) is (/) versus e in the region 
0.0200 < e < 0.0302; (b) is ln(ln(/)25 — ln(/o)25) versus In jet — e| after tangent bifurcation (where 
ln(/o)25 = 2.34- • •); (c) shows the magnified view of the inset in (b); (d) is ln(/)i versus In \et — e\ 
before tangent bifurcation. The average length of the phase locking state is the average rotation 
numbers such that {l)]\f = {T)/2itN, where (T) is the average phase locking time. 



et — e < but also in the region et — e > 0. 



V. DISCUSSIONS 

In the previous study, Lee et a/.^^ once showed that the characteristic relation of the 
average length of the phase locking state for (e^ — e) > (e^ — e)^^/^ that is the scaling of 
type-I intermittency. And then they claimed that the relation deforms to the scaling of eyelet 
intermittency for (e^ — e) < based on the numerical data only in the narrow range 0.0276 < 
e < 0.0286. They obtained the tangent bifurcation point et = 0.0276 and the critical point 
for PS ec = 0.0286 neglecting all the fast fluctuation terms, which are highly important in 
this regard as we have explained above. Also a monograph claimed that the critical point for 
PS is ec ~ 0.028 relying based on the Lyapunov exponent analysis 2^. Unlike their claims, 
however, we have showed that the true tangent bifurcation point is ej = 0.0277 and we have 
obtained the characteristic relation in the wider range 0.0200 < e < 0.0302 where it deforms 
from (/) oc (et — e)^^^"^ to (/) oc exp(a;|et — e|^/^) continuously as e crosses et = 0.0277. This 
deformation is the typical characteristic of type-I intermittency in the presence of noise, as 



3- 
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In numerical simulations, we have 



it was confirmed experimentally in our recent paper 
also found that phase jumps still occur even at e = 0.0304. As mentioned in the above, ^„ 
is bounded in the coupled Rossler oscillators. Nevertheless, it is very hard to determine the 
correct ec because it takes too long time to find out the point where the average length of 
the phase locking state becomes infinite due to the exponential increment of (l). Instead of 
the PS point, the result of numerical simulation follows well Eq. (8) in the region which we 
studied. 

In the route to PS, the phase slip phenomenon in coupled chaotic oscillators is usually 



n coupled cm 

BBQQ. 



described by the Langevin equation: = + C I^M, llil ll^- We understand that the 

■osing the stability for the fixed point (p* where — 4t^L. = by 



phenomenon is in a process of 

the stochastic perturbation ^|, |21|, and so the origin of the transition seems to be univer- 
sal in coupled chaotic oscillators. The characteristic scaling rule can be deformed according 
to the local structure of the Poincare map near the bifurcation point i.e., — Recently 
we observed a similar transition route which has the same origin in coupled hyper-chaotic 
Rossler oscillators whose characteristic scaling rule is governed by type-II intermittency 
in the presence of noise because its normal form has cubic polynomial type instead of a 



quadratic one 



VI. CONCLUSIONS 



We have studied the origin of the transition to PS via phase jumps in coupled Rossler 
oscillators analytically as well as numerically. Analysis of the phase equation and the numer- 
ically constructed local Poincare map reveal that the transition to PS via 2tt phase jumps is 
governed by type-I intermittency in the presence of external additive noise. The character- 
istic behavior of the average length of the phase locking state with respect to the parameter 
e obtained both by numerical fitting and the FPE approach obeys (/) = (/q) exp(Q;|e — |^/^) 
for e > et, with the well known scaling form (Z) oc 1/ a/Ic — | for e < et- 



10 



VII. ACKNOWLEDGEMENTS 

Authors thank S.-Y. Lee for helpful discussions. This work is supported by Creative 
Research Initiatives of the Korean Ministry of Science and Technology. 



[1] L. Fabiny, P. Colet, and R. Roy, Phys. Rev. A 47, 4287 (1993); J.F. Heagy T.L. Caroll and 
L.M. Pecora, Phys. Rev. A 50, 1874 (1994); S.K. Han, C. Kurrer, and Y. Kuramoto, Phys. 
Rev. Lett. 75, 3190 (1995). 

[2] A. Pikovsky, M. Rosenblum, and J. Kurths, Phase Synchronization in Regular and Chaotic 
Systems: a Tutorial (1999). 

[3] S.H. Strogatz, Nonlinear Dynamics and Chaos (Perseus Books, New York, 1994). 

[4] L.M. Pecora and T.L. Carroll, Phys. Rev. Lett. 64, 821 (1990). 

[5] S. Rim, D.U. Hwang, I. Kim, and CM. Kim, Phys. Rev. Lett. 85, 2304 (2000). 

[6] M. Rosenblum, A. Pikovsky and J. Kurths, Phys. Rev. Lett. 76, 1804 (1996). 

[7] A.S. Pikovsky M. Rosenblum, G.V. Osipov, and J. Kurths, Physica D 104, 219 (1997). 

[8] E. Rosa, Jr., E. Ott, and M.H. Hess, Phys. Rev. Lett. 80, 1642 (1998). 

[9] M. Rosenblum, A. Pikovsky and J. Kurths, Phys. Rev. Lett. 78, 4193 (1997). 
[10] K.J. Lee, Y. Kwak, and T.K. Lim, Phys. Rev. Lett. 81, 321 (1998). 
[11] LB. Kim, CM. Kim, W.H. Kye, and Y.J. Park, Phys. Rev. E 62, 8826 (2000). 
[12] A. Pikovsky M. Zaks, M. Rosenblum, G. Osipov, and J. Kurth, Chaos, 7, 680 (1997). 
[13] C. Grebogi, E. Ott, and J.A. Yorke, Phys. Rev. Lett. 50, 935 (1983). 

[14] A. Pikovsky, G. Osipov, M. Rosenblum, M. Zaks, and J. Kurths, Phys. Rev. Lett. 79, 47 
(1997). 

[15] W.H. Kyc and CM. Kim, Phys. Rev. E 62, 6304 (2000). 

[16] J.P. Eckmann, L. Thomas and P. Wittwer, J. Phys. A14, 3153 (1982); J.E. ffirsch, B.A. 

Huberman, and D.J. Scalapino, Phys. Rev. A 25, 519 (1982); J.P. Crutchfield, J.D. Farmer, 

and B.A. Huberman, Phys. Rep. 92, 45 (1982). 
[17] R.L. Stratonovich, Topics in the Theory of Random Noise (Gordon and Breach, New York, 

1963). 

[18] P. Manneville and Y. Pomeau, Phys. Lett. 75 A, 1 (1979). 



11 



[19] J.E. Hirsch, M. Nauenberg, and D.J. Scalapino, Phys. Lett. 87A, 391 (1982); B. Hu and J. 

Rudnick, Phys. Rev. Lett. 48, 1645 (1982). 
[20] CM. Kim, O.J. Kwon, E.K. Lee, and H.Y. Lee, Phys. Rev. Lett. 73, 525 (1994); CM. Kim, 

G.S. Yim, J.W. Ryu, and Y.J. Park, Phys. Rev. Lett. 80, 5317 (1998); O.J. Kwon, CM. Kim, 

E.K. Lee, and H. Lee, Phys. Rev. E 53, 1253 (1996). 
[21] C. W. Gardiner, Handbook of Stochastic Methods (Springer- Verlag, New York, 1985), Second 

Edition. 

[22] H. Risken, The Fokker-Planck Equation (Springer- Verlag, New York, 1996) second edition. 

[23] A. Pikovsky, M. Rosenbhim, and J. Kurth, Synchronization a Universal Concept in Nonlinear 
sciences, (Cambridge University Press, 2001). 

[24] J.H. Cho, M.S. Ko, Y.J. Park, and CM. Kim, Phys. Rev. E 65, 036222 (2002). 

[25] The scahng rule near the tangent point can be estimated as follows. Average laminar length 
is the function of coupling strength such that (l) ~ /(|et — e|). If \et — e| ^ 1 then /(je^ — e|) ~ 
/(O) + rmet - e| + 0{\et - ep). It means that f{\et - e|) « /(O) exp(/'(0)//(0)|et - e|). Thus 
average laminar length is given by (Z) ^ exp(Q!|et — e|) for |et — e| <C 1, where a is a constant. 



